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ABSTRACT 

The time asymptotic behaviour of a relativistic (parallel) 
shock wave significantly modified by the diffusive accelera- 
tion of cosmic-rays is investigated by means of relativistic 
hydrodynamical equations for both the cosmic-rays and thermal 
gas. The form of the shock structure equation and the disper- 
sion relation for both long and short wavelength waves in the 
system are obtained. The dependence of the shock accelera- 
tion efficiency on the upstream fluid speed, long wavelength 
Mach number and the ratio N = P C o/( p co + p go) ( p co and p go are 
the upstream cosmic-ray and thermal gas pressures respectively) 
are studied. 

1. Introduction . Nonlinear theories of diffusive shock acceleration of 
cosmic-rays were initially developed by Axford et al_. (1977, 1982), Drury 
and Volk (1981) (see Drury 1983 for a review and further references). 
These authors discuss how the efficient acceleration of cosmic-rays in 
shock waves by the first order Fermi mechanism leads to shock broadening 
or shock structure due to the interaction between the background thermal 
gas arid cosmic-rays. Previous work on this problem has been restricted 
to non-relativistic shocks. The purpose of the present paper is to ex- 
tend these models to relativistic shocks in which the background fluid 
velocity is an appreciable fraction of the speed of light. 

2. Equations . Consider a one dimensional, steady state inviscid, hydro- 
dynamical model of cosmic-ray acceleration in a (parallel) relativistic 


shock. The basic conservation laws for mass momentum and energy are: 

^ tv V ■= T,v (mass conservation), (1) 

C^i) v i* m (momentum), (2) 

V* V (W c * W 3 ) + (?. Y*- l) t£ p » F*, (energy), (3) 

where J n , P n , F n are the (constant) mass, momentum, and energy fluxes of 
the system, and 

W„=. P c+ - f c - . (4) 

w 3 = * f, - V, P, /6i -j) + vw c*-, (5) 


are the cosmic-ray and gas enthalpy respectively. Pr» are the 

cosmic-ray pressure, energy density and adiabatic index; whereas Pq, 
are the corresponding gas quantities; n denotes the gas density, T«* rep- 
resents the cosmic-rav energy flux in the scattering frame (or fluid frame) 
and Y= (1 - V^/c^)-^ is the Lorentz factor of the fluid speed V. 
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Assuming there are no dissipative mechanisms for the gas in the fluid 
frame, it follows that the thermal gas is compressed adiabatically fol- 
lowing the flow (except at gas subshocks) so that P g = Pg 0 (n/n 0 ) *», 
where the subscript 'o' denotes the (constant) upstream state. To com- 
plete the system we use the cosmic-ray momentum equation (in the scat- 
tering frame): 

-•V -C = of <Afc/^ + K+ « 

governing the interaction between the cosmic-rays and gas via the scat- 
tering frequency V . (Note the hydrodynamical diffusion coefficient in 
this development is K = c.*(Y c -i) / -V. 

3. The Shock Structure Equation . The cosmic-ray momentum equation (6) 
used in conjunction with equations (1) - (5), reduces to the shock 
structure equation 

= c n l (j'r *50} (?) 

governing the inverse compression ratio y = V V/y^of the flow. 

In (7), can be expressed solely in terms of y (or V) using (2) and 
(3), and a = (y f /W)' 1 *- denotes the nonrelativistic thermal gas sound 
speed. 


Note in particular the singularity in the denominator in (7) occurs 
when the fluid speed matches the short wavelength thermal gas sound speed 

V s =■ ^ ^ • (8) 

For non-relativistic shocks ; also note that there is a strict up- 

per limit to the thermal gas sound speed V s in the relativistic case with 
V s <CViTc.0n the other hand, for a near uniform shock transition we 
have T * 0 throughout the flow; the solution of T*° = 0 for the up- 
stream fluid speed V 0 then yields the long wavelength sound speed in the 
system 

X . = +“>•)/ 1 0 ) 

where a c ', = ?*./«.„ . 


A further complication that arises in the relativistic theory (as op- 
posed to the non-relativistic theory) is that T^° in (7) develops a sin- 
gularity when the fluid speed t which is the sound speed of the 

‘effectively massless' cosmic-ray gas (cf. Weinberg, 1972). The occur- 
rence of this singularity is presumably related to the fact that the adi- 


abatic deceleration rate for a fully relativistic cosmic-ray gas is given 
by 

where ^ and p 1 and m‘ represent the individual particle 


momentum and relativistic mass in the scattering frame. The basic impli- 
cation of (10) i s th at cosmic-rays are accelerated in compressive flows 
with V < s 4V C -T c 5 U t are decelerated in compressive flows with 
v> c/\Ti (cf. Webb, 1985). 
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4. Numerical Results. The nature of the shock transition can be ana- 
lyzed by a geometrical method employed by Drury and Volk (1981) in thei 
work on non-relativistic hydrodynamical models. This is illustrated in 
Figure 1 for a shock in which the long wavelength Mach number M = 2 
N = P co /(Pco + Pgo) = 0.3, ; z 0 = V 0 /c = 0.5, ^ = 4/3 and X = 5/3 
shows the Hugoniot, gas adiabat, short wavelength some line and 
curve in the (P q , V/c) plane. The Hugoniot corresponds to the locus 
T'** = 0 in the y (P q , V/c) plane, whereas the short wavelength sonic line 
corresponds to the zero of the denominator in (7). The upstream state 
(A) and downstream state (C) lie on the Hugoniot. In the initial part 
of the transition the gas is compressed adiabatically (AB) followed by a 
subshock to a uniform downstream state (the straight line segment BC). 
Across the subshock the momentum and energy fluxes of the cosmic-rays 
and thermal gas are balanced separately. Figure 2 shows the variation 
of the cosmic-ray shock acceleration efficiency *?, as a function of the 
upstream fluid velocity Vo for a range of upstream Mach numbers (M = 2, 

4. 5, 10), N = 0.3, K c = 4/3 and Kg = 5/3 ( is defined as the fraction 

of flow kinetic energy that is converted to cosmic-ray, energy in the 
transition, i.e., vj t /\f li l VVA/ 0.5773, there is 

a weak dependence of ») s on z 0 , and vj t increases monotonically with M. At 
0.577 < z 0 <0.9, there are no solutions and the solutions with z 0 > 0.9 
are of low efficiency. The negative efficiency solutions at z 0 > 0.9 are 
obtained in expanding flows (V increases ini tially) followed by a sub- 
shock to a downstream state z 2 = V 2 / c ^ v V c -| . In compressive flows with 
V 0 > c (^ fc - l) 1 ' 2 = 0.577c, the cosmic-ray pressure decreases as the thermal 
gas is compressed adiabatically, and the cosmic-ray pressure may even- 
tually become negative; it is then necessary to insert a subshock in the 
flow while P c > 0; this behaviour is in accord with the result (10). 

Figure 3 shows the variation of the downstream cosmic-ray pressure P C 2 as 
a function of the downstream cosmic-ray adiabatic index ^c % ( \ 0 = 1-55, » ^ — 
5/3, M = 5, z 0 = 0.5) for a range of N. The dotted curves correspond to 
smooth transition solutions; whereas solutions with subshocks are given by 
the full (solid) curves (cf.Achterberg et al_. 1984; Heavens 1984 for the 
non-relativisitc case). 

5. Conclusions . The main point to emerge from the present study is that 
to efficiently accelerate cosmic-rays in relativis tic shocks, the up- 
stream fluid speed V 0 needs to be less than V cr = jY c -l c; compressive 
flows with V>V cr lead to a loss of internal energy to the C.R. gas. 
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